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ABSTRACT 

We present a preliminary survey of 58 radio sources within the isoplanatic patches (r < 25") of 
bright (11 < R < 12) stars suitable for use as natural guide stars with high-order adaptive optics 
(AO). An optical and near-infrared imaging survey was conducted utilizing tip-tilt corrections in the 
optical and AO in the near-infrared. Spectral Energy Distributions (SEDs) were fit to the multi- 
band data for the purpose of obtaining photometric redshifts using the Hyperz code ( jBolzonella et al. 
2000|. Several of these photometric redshifts were confirmed with spectroscopy, a result that gives 
more confidence to the redshift distribution for the whole sample. Additional long- wavelength data 
from Spitzer, SCUBA, SHARC2, and VLA supplement the optical and near-infrared data. We find 
the sample generally follows and extends the magnitude-redshift relation found for more powerful 
local radio galaxies. The survey has identified several reasonably bright (H = 19 — 20) objects at 
significant redshifts (z > 1) that are now within the capabilities of the current generation of AO-fed 
integral-field spectrographs. These objects constitute a unique sample that can be used for detailed 
ground-based AO studies of galactic structure, evolution, and AGN formation at high redshift. 
Subject headings: galaxies: high-redshift — galaxies: fundamental parameters — galaxies: evolution 
— galaxies: formation 



1. INTRODUCTION 



One of the most active debates in extragalactic as- 
tronomy is on the nature of spheroid formation at high 
redshifts. The hierarchical galaxy formation scenario is 
base d on the conclusions of semi-analytic ACDM mod- 



els (|Thomas fc Ka uffmann 1999] |Kauffmann et al ||1999 
Kauffmann fc Haehnelt 2000) as well as observations oi' 



galaxies (Miyazaki et al. 2003 Dickinson et al.| 
van Dokkum et al. 2008)/ Itis typically describe 
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J . It~is typically described as 
a building up of a large stellar system from smaller ones, 
and predicts that massive early-type galaxies should 
have undergone final assembly at relatively late epochs 
(1 < z < 2). Conversely, the monolithic forma tion sce- 
nario models, supported by other observations dThomas 
et al.|2002||Brodie fc Strader|2006|[McGrath et al.|2008p , 
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require that massive spheroids form at much higher red- 
shifts directly from primordial density fluctuations (see 



Peebles (2002) for a review and discussion) and evolve 
passively to the present epoch. 

In order to test the predictions of these scenarios, 
and thereby determine which provides the more accu- 
rate model of galaxy formation, we have attempted to 
study massive galaxies at high redshift by observing the 
host galaxies of radio sources. Though these objects are 
not classified as normal galaxies, this technique has cer- 
tain advantages over other high redshift galaxy search 
methods (e.g. Lyman-break, and NIR-selection). High 
redshift radio galaxies (HzRGs) tend to have large stellar 
masses and high near-infrared luminosities making them 
accessible to ground-based observations. Furthermore, 
radio galaxies follow a magnitude -redshift relat i on (th e 
near-IR Hubble diagram, see Ivan Breugel et all (1 19981)), 
which allows a selection of high redshift objects based 
on their apparent brightnesses. Lilly & Longair (1984]) 
also found that the K-band magnitude-redshift (K-z) re- 
lation of the 3CR catalog of powerful radio galaxies could 
be well fit by a passively evolving old stellar population, 
similar to present-day elliptical galaxies. The tendency 
for HzRGs to reside in over-dense regions in the early 
universe, makes them likely progenito rs of the brighte st 



cluster galaxies in the present epoch (Bes t et al.|[l998| ). 
There are a few disadvant ages to observing HzRGs. 



First, the "alignment effect" (|Chambers et al.| 1 987; Mc- 
Carthy et al.||1987| [Chambers fc Miley||1990| |Pentericci 
et al. 200 ip was identified when it was discovered that 



2 



these objects' optical morphologies at high redshift were 
aligned with structures in the radio, and therefore their 
optical/infrared magnitudes could not be treated as in- 
dependent of their radio properties. Several mechanisms 
have been proposed to explain the alignment effect (e.g., 
anisotropic interactions between cluster members, jet- 
triggered star formation, optical synchrotron emission, 
inverse Compton scattering of the CMB photons, and 
thermal conti nuum emiss i on fro m the plasma ionized by 
the AGN, see |McCarthy| fll993| )). Each of these mecha- 
nisms has been shown to be present in certain objects, 
but because no single explanation is satisfactory for all 
cases, it is generally thought that two or more may be the 
source of the alignment effect. Regardless of the causes 
of the alignment effect, in order to obtain the fundamen- 
tal parameters of the underlying host galaxy, the con- 
tributions from the AGN and the host galaxy must be 
disentangled. 

Studies of the relative contribution of evolved stellar 
population and flat-spectrum (and presumably aligned) 
components in optical and infrared images of 3CR ra- 
dio galaxies revealed that the components responsible 
for the alignment effect contribut e only about 10% to 
the total SED of the host galaxy QBest et al.||1998[ ) af- 
fecting both continuum and emission-line morphologies 
( |McCarthy et al. 1987). The alignment effect also gen- 
erally diminishes at wave lengths longward of the 4000 A 



break (Rigler et al. 1992), although there are some ex- 
ceptio ns (Eisenhamt &; Chokshi||l990} | Chambers et al. 



19881). When the alignment effect is examined in de- 
tail, only one-fourth to one-third of low redshift radio 
galaxies have any detectable morphological peculiarities 
(at the fly > 25mag/arcsecond 2 level), and this fr action 



becomes smaller in less powerful radio ga laxies ( Heck- 
|man et aT||1986| |Dunlop fc Peacock| [l993), as a weaker 
radio source has less of an impact on the rest of the 
galaxy. This suggests the millijansky-level radio source 
population predominantly consists of sources with radio 
fluxes sufficiently low that the optical/near-IR morpholo- 
gies and SEDs of host galaxies are probably completely 
dominated by the stellar population of the host galaxy. 

An additional disadvantage of observing high redshift 
galaxies results from the loss of angular and spatial 
resolution when using traditional ground-based obser- 
vational techniques, which is needed to accurately de- 
termine the galaxy's fundamental properties. This can 
easily be achieved from space using HST, though in 
the near-infrared (where these objects' SEDs are less- 
contaminated by their AGNs) HST is limited by both 
the diffraction size (about 0.2" in K-band) and light- 
collecting aperture. This makes a survey of these sources 
extremely time-consuming on a highly competitive tele- 
scope and spectroscopy on these faint sources is almost 
impossible. However, with the advent of ground-based 
AO techniques, a large imaging survey can be conducted 
with a comparable spatial resolution and exposure time 
on a 3-4 meter class telescope, with the limitation that 
most AO systems require a bright guide star in proximity 
to the galaxy on the sky. Luckily radio sources are com- 
mon enough that a significant number meet this criteria 
in wide-area radio survey data sets. 
The Faint Images of the Radio Sky at T wenty- 



centimeters (FIRST) Survey (White et al. 1997) at the 



Very Large Array (VLA) has 80% completeness down 
to 1 mJy and better than \" astrometry that can be 
used to search for HzRGs. The catalog is only mini- 
mal l y cont aminated (about 10% by number, see |Jack-| 
|son| ( |2005| )) by low redshift star forming galaxies which 
quickly become more numerous at levels fainter than 0.1 
mJy. More importantly, the VLA's astrometry allows 
for mostly straightforward optical counterpart identifica- 
tions. And lastly, the density of sources (about 90 sources 
per square degree over the survey's 10,000 square de- 
grees) is sufficiently high that a sizable sample of sources 
are within the isoplanatic patch (about 25" in the NIR) 
of a V<12 star. With these criteria, this survey pro- 
vides a highly efficient means to preselect likely candi- 
dates for high redshift galaxies with undisturbed stellar 
populations and light profiles that can be observed at 
high spatial resolution from the ground initially using 
natural guide star AO. As laser guide star and multi- 
conjugate AO (which also require nearby natural guide 
stars) become more mature, this sample would also be 
well-suited for these instruments and as an initial target 
list for JWST as it relates to several of its key science 
objectives. 

In this paper, we present the first phase of a project, 
involving an imaging survey plus a supplemental spectro- 
scopic survey of the FIRST-BNGS sample. The imaging 
survey consists of optical, NIR, FIR, submillimeter, and 
radio observations. However, the focus of the imaging 
survey was on the optical and NIR wavelength regions 
(since we are mainly interested in the stellar populations 
of these objects). The goal of this survey was to obtain 
multi-wavelength photometry to provide identifications 
and probable redshifts for candidates for high precision 
diagnostics of galaxies at high redshift to study galactic 
structure, evolution, and AGN formation at high red- 
shift. Supplemental spectroscopic observations were also 
obtained to refine redshifts or remove ambiguities in the 
photometric redshift solutions. In section 2, we intro- 
duce the sample we compiled to search for high redshift 
radio galaxies. Section 3 describes our observations and 
reduction process for the imaging survey. In section 4, 
we present our photometric measurement and redshift 
fitting procedures and results. Section 5 discusses the 
spectroscopic observations, reduction process, and re- 
sults, while section 6 briefly summarizes the VLA radio 
data. Finally, in section 7, we summarize our findings 
and discuss the potential of this sample. 

For this paper we adopt (unless otherwise stated) 
the Friedmann-Lemaitre world cosmological model with 
^o=0.3, ^a=0.7, and #o=70km/s/Mpc, giving the 
present age of the Universe as 13.47 Gyr. Also, the 
spectral index, a, will be defined such that the flux of 
a source, S U1 is proportional to v a . 

2. SAMPLE SELECTION 

A cross-correlation of the VLA FIRST surv ey and the 
USNO-A2.0 Catalog ( [Monet fc et^[l998l ) yields 58 
sources with S\agh z > 1-mJy, galactic latitude \bu\ > 35 
located within an annulus 15 < r < 25 arcseconds 
around a 11 < R < 12 star. These criteria ensure that 
the FIRST source lies within the isoplanatic patch of a 
sufficiently bright guide star for the NIR AO observa- 
tions under most seeing conditions without introducing 
additional observing or photometric complications (e.g., 
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PSF wings of the guide star contaminating the sky back- 
ground if the source is too close to the star or the core of 
the guide star saturating significant sections of the detec- 
tors if the star is too bright). We expect n ear diffraction- 
limited performance longward of 1.2/ira ( |Roddier 1999) 
for a typical high-order AO system at Mauna Keal m 
the optical, where only tip-tilt correction is available, the 
point-source sensitivity is increased by up to a factor of 
four over non-compensated images. Since we are proba- 
bly not resolving the high redshift (z>l) FIRST sources, 
we would expect similar performance. 

Table 1 lists the coordinates of each radio source and 
corresponding USNO star corrected for proper motion 
to 2005.0. These 58 objects comprise the FIRST-BNGS 
sample used in this paper . About 90% are expected to 
be FRI or FRII galaxies ( |Jackson||2005| ), and based on 
scaling the 151MHz local radio luminosity function to 1.4 
GHz, probably half will be at significant redshift (z > 1). 

3. OBSERVATIONS 

Here we describe our optical and near-infrared imaging 
data. Figures 1-8 show 54"x54" FIRST radio contour 
maps and 14 // xl4 // thumbnail NIR images for each source 
overlaid with FIRST radio contours. 

3.1. Optical Imaging 

The broadband optical imaging data of the FIRST- 
BNGS sample were obtained with the Or thogonal Par- 
allel Transfer Imaging Camera (OPTIC; |Tonry et al. 
(2002)) mounted at the f/10 focus of the University of 
Hawari 2.2-meter (UH2.2m) telescope and at the Nas- 
myth focus of the Wisconsin, Indiana, Yale, NOAO 
(WIYN) 3-meter telescope at Kitt Peak National Obser- 
vatory. This camera utilizes an effective "tip-tilt" correc- 
tion feature of a specialized CCD, called an orthogonal 
transfer array, with no physically moving parts. This 
technique moves the accumulated charge from an astro- 
nomical source around on the CCD based on the centroid 
of a guide star read out at short intervals (20-100ms) in 
a nearby region of the chip. All of our sources are well 
within a BNGS isoplanatic patch for tip-tilt correction 
(about 10' for a 2 meter telescope). The tip-tilt usually 
improves the point-source sensitivity by nearly a factor of 
four in most seeing conditions and slightly increases the 
achieved spatial resolution. Our best measured FWHM 
was 0.3" for a 300 second exposure. Since we expect the 
angular size of these high redshift sources to be about 
this scale, tip-tilt correction only provides an increase in 
efficiency for flux measurements. 

Our observations consisted of several sequential or- 
thogonally transferred (OT) 300-second exposures in var- 
ious broadband filters (B, V, R, I, and z'). Data were 
taken on about 50 nights between 2002 December and 
2005 April (Table 2) under photometric conditions. The 
6' field of view was oriented strategically with the bright 
guide star in the guiding region of the chip and the object 
in the imaging region. Due to the proximity of the sci- 
ence target to the guide star, diffraction spikes and PSF 
wings were sometimes a problem. In addition, those ob- 
jects that were particularly close to the guide star were at 
the edge of the detector's science region. As the guiding 
regions are at the top and bottom of the detector, the in- 
strument was occasionally rotated 90 degrees in order to 
accommodate target objects to the east and west of the 



guide star. Supplemental rough guiding was also accom- 
plished by OPTIC communicating frequent guide offsets 
to the telescope control system. 

With the exception of the generation of flat-field frames 
("flats"), standard IRAF procedures were for the data re- 
duction. The flats were constructed using a special pro- 
gram, kindly provided by John Tonry, called "conflat", 
which constructs a flat field from any normal flat image 
(in this case, the flat was created from the median of sev- 
eral high-signal dome flat images) by convolving it with 
the sequence of OT shifts executed by OPTIC during 
the science exposure. A separate convolved flat is there- 
fore needed for each exposure. Each set of images was 
flat-fielded, background subtracted, aligned using bright 
stars, then averaged together. An absolute position was 
derived from the locatio n of the bright guide s tar using 
the USNO-A2.0 catalog ([Monet fc et al.||1998[ ). A sim- 
ilar procedure was used f or sequences ol short exposure 
on photometric standard Landolt (1992) fields, but the 
frames were not stacked so that the individual exposures 
could be used to estimate uncertainties in the calibration. 
Photometric zeropoints were derived based on these Lan- 
dolt standards in the Vega magnitude system. 

3.2. NIR Imaging 

For the NIR survey of 58 FIRST-BNGS sources, our 
strategy was to obtain at least H-band photometry for 
the entire sample. Because the HzRGs in the sample 
are relatively compact, this was most efficiently accom- 
plished using the 3.6-meter Canada- France- Hawai'i Tele- 
scope (CFHT) wit h the the Pueo AO b onnette and KIR 
infrared detector (Rigaut et al. 1998). Pueo-KIR in- 



corpo rates a wavefront curvature sensor (Roddier et al. 
1991) and a 19-electrode deformable mirror! This pro- 
vides near diffraction-limited imaging in good conditions 
using guide stars similar to those in the FIRST-BNGS 
sample. Since the diffraction limit of CFHT is around 
0.2" in H-band, most the high redshift sources are barely 
resolved or unresolved, and therefore we are presenting 
only flux measurement from these observations. How- 
ever, the increase in point-source sensitivity over non- 
AO NIR observations allowed us to rapidly observe all 
sources in the sample in a handful of nights. Some addi- 
tional non-AO J, H, and K photometric d ata were taken 
with the Quick Infrared Camera (QUIRC; Hodapp et al. 
( 1995 ])) on the UH 2. 2m telescope and the SpeX imager 
(Rayner et al. 2003) on the IRTF to supplement the 
CFHT H-band data. A subsequent K-band AO imaging 
campaign to measure morphologies of these sources was 
carried out on a subsample of 18 FIRST-BNGS high red- 
shift candidates (Stalder & Chambers, in prep.). These 
observations were made using the Subaru 8-meter tele- 
scope and the Infrared C amera and Spectrograph (IRCS) 
( Kobayashi et al.||2000 ) mounted behind the 36-element 
curvature- sensing AO system and the photometry is in- 
cluded in this data set. 

The J and H AO data were collected under excel- 
lent seeing and photometric conditions over the course 
of 19 nights (Table 2). The K-band AO data were ob- 
tained during photometric conditions, but seeing con- 
ditions were below average (1.5-2.0 // natural seeing in 
V-band) on several of the nights, unfortunately giving 
a mediocre corrected K-band FWHM (0.2-0.5"). The 
non-AO data were under photometric and varying see- 
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Object Positions (J2000.0). 
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FlG. 1. — Images of 8 radio source fields from the BNGS sample. Images are centered on the optical identification or center of the radio source if 
undetected. The left panel is the FIRST contour radio map. The right panel is H-band Pueo-KIR AO images. The FIRST radio map is 54" on a 
side (the contours are set as 0.25, 0.5, 1, 2, 4, 8, 16, and 32 mJy) and the H-band images are 14 7/ on a side. The contours are set to be 5 equal linear 
levels from the peak flux in the field to mJy flux density. At left, from top to bottom shows the optical/IR sources, F0023-0904, F0129-0140, 
F0152-0029, and F0152+0052. At right, from top to bottom shows, F0202-0021, F0216+0038, F0916+1134, and F0919+1007. Note that there is 
the contamination from negative (white) amplifier and positive (black hourglass) ghosting artifacts (in Pueo-KIR images) from the saturated guide 
star. The edges of the detectors are also visible in some frames. 




FlG. 2. — Same as Figure 1, for 8 additional objects. This set at left, from top to bottom shows, F0938+2326, F0939-0128, F0942+1520, and 
F0943-0327. This set at right, from top to bottom shows, F0950+1619, F0952+2405, F0955+2951, and F0955+0113. 



ing conditions. 

Calculated Strehl ratios (the ratio of the maximum 
of the on-axis PSF to the maximum of a theoretical 
diffraction-limited PSF) for Pueo-KIR were around 10% 
in J and consistently above 40% and reached as high as 
80%-90% in H, and remained stable throughout these 
nights. The IRCS K-band Strehls were around 20-25%. 
These are typical Subaru Strehl ratios because the sys- 
tem has only 36 subapertures on an 8-meter diameter, 
whereas Pueo has 19 subapertures over the 3.6-meter di- 



ameter that more finely samples the isoplanatic regions 
over the aperture. Both AO systems resulted in near 
diffraction-limited performance in average to good nat- 
ural seeing (0.15" on-axis corrected FWHM for CFHT 
in H and 0.08" for Subaru in K). Unfortunately, due to 
worse than average seeing conditions, little of the IRCS 
data achieved this. The exposure times were 180-seconds 
in J, 120-seconds in H, and 90 in K. A custom dithering 
pattern was used that kept the bright guide star on the 
field for aligning the images, but did not allow the regions 
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FlG. 3. — Same as Figure 1, for 8 additional objects. This set at left, from top to bottom shows, F0956-0533, F0958+2721, F1000-0636, and 
F1008-0605. This set at right, from top to bottom shows, F1010+2527, F1010+2727, F1014+1438, and F1016+1513. 



saturated by the guide star to overlap with the object, 
to avoid persistence problems in photometric measure- 
ments and in flat frame construction. For the Subaru 
AO data, an off- axis PSF star at the same distance to 
the guide star and similar brightness to our science tar- 
gets was observed each night to represent the PSF in our 
fitting process for the K-band data. 

The Pueo-KIR detector suffers from additional ampli- 
fier related artifacts from the saturated guide star re- 
flected in each quadrant, plus a "ghost" guide star due 



to internal reflections within the instrument. These were 
easily identified because of their consistency, but are 
remarked on here because of the effect on the overall 
quality of the data. The amplifier artifacts can be seen 
as negative (white) signals in most H-band thumbnail 
images, for example, just southwest of the F0202-0021 
object at H-band in Figure 2. The "ghost" image can 
plainly be seen as an hourglass shaped object to the east 
of F02 16+0038 in Figure 2 and also in the H-band image 
of F1010+2727 in Figure 3. Some objects were also near 
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FlG. 4. — Same as Figure 1, for 8 additional objects with the exception that F1140+1316 is OPTIC I-band with the same spatial scale and 
FIRST contours as the H-band Pueo images. This set at left, from top to bottom shows, F1024-0031, F1027+0520, F1039+2602, and F1040+2323. 
This set at right, from top to bottom shows, F1116+0235, F1133+0312, F1140+1316, and F1147+2647. 



the edge of the detector FOV, as seen in several images 
(a consequence of keeping both guide star and science 
object on the chip). 

The non-AO data was taken over 5 nights in 2003 Au- 
gust using QUIRC on the UH2.2m and 2006 August 
using SpeX on IRTF with similar exposure times and 
dithering patterns as the CFHT Pueo-KIR and Subaru 
IRCS data. 

Standard IRAF procedures were used to reduce the 



near-infrared data. Median sky flats were constructed 
from all of the exposures of each field then normalized 
by the mode of the sky values. For each exposure, the 
data was divided by the normalized flat and sky sub- 
tracted based on the mode. Finally, the science frames 
were stacked, using the bright guide star for registration. 
World coordinate system data was also entered into the 
file header based on the position of the guide star. The 
same procedure was used for short exposure photometric 
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FlG. 5. — Same as Figure 1, for 8 additional objects. This set at left, from top to bottom shows, F1155+2620, F1158+1716, F1202+0654, and 
F1211+3616. This set at right, from top to bottom shows, F1215+3242, F1217-0529, F1217+3810, and F1218-0625. 



standard fields (by using the science object's median sky 
frame closest in time to when the sequence was taken), 
but the fields were not stacked in order to obtain multi- 
ple measurements. Photometric zeropoints were derived 
from the magnitudes in th e UKIRT Faint Standard cat- 
alog ( Leggett et al. 2006) for the Mauna Kea filter set 
(|Sinions & Tokun aga 2002} |Tokunaga et al.|2002| ) on the 
Vega magnitude system. 

3.3. Spitzer Archive Data 



A search of all 58 FIRST-BNGS sources in the Spitzer 
Archives resulted in serendipitous observations of 2 ob- 
jects (F1237+1141 and F1430+3557). The Spitzer IRAC 
3.6/5.8/im and MIPS 24/im, 70/im, and 160/im mosaiced 
images were produced from basic calibrated data (BCD) 
images using the Mosaick ing and Point S ource Extrac- 
tion tool, MOPEX QMakovoz fc Khan|2005| . Some clean- 
ing was done for the MIPS 70 and 160/im data also using 
MOPEX. Aperture photometry was accomplished using 
the IRAF "phot" routines with an appropriate aperture 
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FlG. 6. — Same as Figure 1, for 8 additional objects with the exception that F1329+1748 is OPTIC I-band with the same spatial scale and 
FIRST contours as the H-band Pueo images. This set at left, from top to bottom shows, F1218-0716, F1234+2001, F1237+1141, and F1315+4438. 
This set at right, from top to bottom shows, F1329+1748, F1355+3607, F1430+3557, and F1435-0019. 



correction from the respective IRAC and MIPS websites. 
Sensitivity was measured through fitting the background 
assuming Poisson statistics. For F1237+1141 at 3.6/im, 
the guide star has some diffraction spikes near the object, 
but was mitigated by measuring the flux with a smaller 
aperture and separately at two observed rotation angles 
(with different spike orientations). For 5.8/im and the 
MIPS channels, the nearby guide star was much dimmer 
at these wavelengths so the flux measurements were not 



affected by gradients or artifacts. 

3.4. Submillimeter Bolometer Observations 

In an attempt to detect submillimeter emission from 
the objects in the sample, we obtained time on CSO and 
JCMT. We were able to observe 2 objects in the sam- 
ple using the CSO bolometer array, SHARC2 over the 
course of several nights. SHARC2 data were reduced us- 
ing the CRUSH pipeline (with faint and compact flags 
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FlG. 7. — Same as Figure 1, for 8 additional objects with the exception that F1445+2702 is OPTIC Fband with the same spatial scale and FIRST 
contours as the H-band Pueo images. This set at left, from top to bottom shows, F1445+2702, F1447+1217, F1451+0556, and F1458+4319. This 
set at right, from top to bottom shows, F1505+4457, F1524+5122, F1644+2554, and F2217-0837. 



enabled), and aperture photometry was measured with 
IRAF routines. Though no detections were made, the 
background noise was measured to get sensitivity for each 
field. Several hours of JCMT SCUBA queue time was 
also awarded, though only a few observations were actu- 
ally executed on 2 sources. The photometric data were 
reduced and extracted using the ORAC-DR pipeline us- 
ing the nearest calibrators. A weak detection was made 
of Fl 116+0235, which is a low-redshift source. Table 3 



shows the measurements from these data. 

3.5. Identification 

The reduced and registered optical and near-IR im- 
ages were compared to the FIRST survey cutouts at the 
coordinates of the radio source in the catalog. Optical 
identifications of single radio sources were straightfor- 
ward. We required the counterpart to be within 3" of 
the peak of the radio source. This allowed for the FIRST 
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FlG. 8. — Same as Figure 1, for 2 additional objects. At left is F2217-0138 and at right is F2354-0055. 



TABLE 2 

OPTIC AND NIR OBSERVATIONS 



Telescope 


Instrument 


Filters 


Dates (UT) 


FWHM a 


UH 2.2-meter 


OPTIC 


B,V,R,I,z' 


2002 


April 23-27 


1.0"-1.5" 








2002 


December 7-10 


0.7"-1.2" 








2003 


January 24-27 


0.5"-0.9" 








2003 


March 20-25 


0.6"-1.2" 








2003 


May 20 


1.0"-1.2" 








2003 


August 1-5 


0.7"-1.0" 








2004 


February 27 


0.8 // -1.0 // 








2004 


July 28 


0.5" 








2004 


August 21-24 


0.8 // -2.0 // 








2005 


February 10-17 


0.7 // -2.5 // 


WIYN 


OPTIC 


B,V,R,I,z' 


2003 


October 14 


0.6"-1.2" 








2004 


June 4 


1.0"-1.2" 








2004 


November 9-10 


1.2"-1.5" 








2005 


April 26-28 


1.0 // -1.2 // 








2006 


January 1 


1.0" 








2006 


June 29-31 


0.5 // -2.0 // 


UH 2.2-meter 


QUIRC 


K' 


2003 


August 10-11 


1.0" 


IRTF 


SPEX 


J,H,K 


2006 


August 14-15 


0.7"-0.9" 


CFHT 


Pueo-KIR 


J,H 


2003 


March 14-16 


0.5"-0.7" 








2003 


October 11-14 


0.2"-0.3" 








2004 


April 3-4 


0.3"-0.4" 








2004 


October 1-3 


0.2"-0.4" 








2005 


April 17-23 


0.2"-0.3" 


Subaru 


IRCS 


K' 


2004 


February 2-3 


0.2"-0.5" 








2004 


November 29 


0.5" 








2005 


January 16 


0.4" 








2005 


February 18-19 


0.3"-0.6" 



a Corrected FWHM range in observed band. 
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TABLE 3 

Spitzer, SHARC2 and SCUBA Fluxes and 3-a upper-limits. 



Name 


instrument 


dates(UT) 


A(jiim) 


F,(mJy) 


ov(mjy) 


F0152-0029 


SHARC2 


Novl6-19,2005 


350 


<24.6 




F1116+0235 


SCUBA 


Nov22,2003 


450 


<39.0 




F1116+0235 


SCUBA 


Nov22,2003 


850 


4.1 


2.1 


F1237+1141 


IRAC 


May27,2004&Jun09,2004 


3.6 


0.0673 


0.0013 


F1237+1141 


IRAC 


May27,2004&Jun09,2004 


5.8 


0.1928 


0.0017 


F1237+1141 


MIPS 


Jun22,2004&;Jun23,2004 


24.0 


0.725 


0.073 


F1237+1141 


MIPS 


Jun22,2004&Jun23,2004 


70.0 


3.21 


0.64 


F1237+1141 


MIPS 


Jun22,2004&Jun23,2004 


160.0 


12.5 


2.5 


F1430+3557 


MIPS 


Feb01,2004 


24.0 


0.374 


0.037 


F1430+3557 


MIPS 


Feb01,2004 


160.0 


<33 




F1644+2554 


SCUBA 


Aug27,2003 


450 


<90 




F1644+2554 


SCUBA 


Aug27,2003 


850 


<6.6 




F2217-0138 


SHARC2 


Novl9,2005 


350 


<195 
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positional accuracy (about 1") and for slightly elongated 
radio structures where the central source positions are 
not well measured. For double-lobed radio sources to be 
matched, detections had to be along the radio axis and 
near the center of the two lobes. For irregular radio mor- 
phologies, many objects in the corresponding area were 
marked as possible matches (i.e., N, S, E, W, etc.). The 
objects nearest to the radio source peak were usually se- 
lected as the most probable counterparts. 

3.6. Photometry 

Calibrated flux measurements were made using the x 2 - 
minimization surf ace-brightness fitting routine GALFIT 
( Peng et al.|2Q Q2) for each object for all optical and NIR 
bands observed. GALFIT was chosen due to the den- 
sity of sources and presence of artifacts in the data as 
it fits multiple objects simultaneously and is better able 
to deal with artifacts than aperture photometry. For the 
B, V, R, I, z', J, H and non-AO K no PSF-convolution 
was used (the Pueo-KIR J and H data have relatively 
low S/N and do not have an observed off-axis PSF star), 
so each object was fitted to a pure Sersic (J ~ r 1 / 71 ) pro- 
file which spans the possible profiles of unresolved Gaus- 
sian (n = 0.5), exponential disk (n = 1.0), and de Vau- 
couleurs (n = 4.0) . For the K-band Subaru AO data, 
the objects were fit as Sersic profiles convolved with a 
normalized PSF derived from an off- axis 17th magnitude 
star imaged the same night and if necessary, smoothed to 
the FWHM measured in the science frame. These stars 
were selected to be at the same distance to its guide star 
as the science objects, and at a similar airmass, so should 
capture the long temporal scale off-axis PSF as long as 
conditions did not change much throughout the night. 
This was verified with other point sources in our science 
fields. Because the background is of critical importance 
in these fits, tests were run using different background fit- 
ting algorithms, using both aperture fitting and using the 
GALFIT integrated algorithm. The aperture-measured 
background was found to be more accurate a nd repro- 
ducib l e; the same conclusion was reached by Haussler 
using both real and simulated data. Resid- 
uals were checked to verify that the final fit was good. 
Only total magnitudes were used from these fits. The 
magnitude errors reported by GALFIT were larger than 
the 0.03 mag intrinsi c scatter for S/N > 2.5 found by 
Haussler et al.| ( |2007| ), so the GALFIT photometric er- 
rors added in quadrature to this was used as an estimate 
of the error. The result for each object was total flux 
in each observed filter. These measurements are aper- 
ture and atmospheric-seeing independent as they are de- 
rived from a model fit. This was confirmed by smooth- 
ing sample data and GALFIT recovering essentially the 
same flux (within measurement errors). These photome- 
try measurements were also robust against gradients and 
artifacts caused by the nearby bright guide star, this was 
verified by repeatability of the measurements of the same 
object at various position angles and with the guide star 
on and off the detector. In cases where an artifact is 
particularly close to the science object, they are mod- 
eled in the GALFIT fit along with the science object, 
producing excellent residuals. Table 4 shows the best-fit 
magnitudes from GALFIT or 3-sigma upper-limits based 
on background noise in the data and assuming Poisson 
statistics. 



4. SED FITTING 

4.1. Optical-NIR Photometric Redshifts 

The initial redshift determinations for each object with 
4 or more bands of photometric data (plus 1 object with a 
large break between 2 bands, F0942+1520) w ere derived 
using the p ublic SED fitting co de, Hyperz (|Bolzonella| 
et al.[|2000| . The library of t he |Mannucci et al | j20UTf 
near-infrared extensions to the Kinney et al. (1996) tem- 
plates for E, SO, Sa, Sb, Sc, and SB1 to SB6 (starburst 
models over a range of reddening) were used as the range 
of possible SEDs for this sample of radio sources (see Ap- 
pendix A for a more detailed description and comparison 
plots of the templates). Extensions in the UV and mid 
to far IR wavelength regions were also added (see Ap- 
pendix). 

Hyperz uses a ^-minimization procedure to determine 
the best-fit SED and redshift after applying a c orrection 
for Lyman forest absorption from Madau (1995). Within 
the limitations of the code, the redshift range was set to 
between < z < 7 in steps of 0.05, and the extinction 
in the rest-frame V-band {Ay) was allowed to vary from 
to 10 magnitudes. Reddening was calculated using the 



Calzetti et al. pOOO ) extinction law. Figures 9-17 shows 
the best-fit SEDs for these objects with the parameters 
listed in Table 5. 

Since only low redshift templates were used, redden- 
ing, which is fit for by Hyperz, was relied on to find 
most likely solution for a given object at each redshift. 
This is reasonable given the degeneracies between star 
formation history, age, metallicity and reddening. It 
should be pointed out that therefore no detailed infor- 
mation beyond redshift, absolute magnitude, and general 
SED-type (early, late, star-burst) about these parame- 
ters should be considered certain, so were not given in 
our analysis. 

All fits were confirmed using our own software that car- 
ries out a procedure similar to that of Hyperz. The code 
performs a x 2 - mm i m i za tion to fit the SED which first 
converts magnitudes to flux, then creates a library of ev- 
ery SED template at ev ery redshift ste p in the given red- 
shift range, applies the |Madau| (1995 ) absorption curve 
to them, integrates each redshift ed model over the filter 
bandpasses, computes the least squares fit between the 
observed fluxes and the model fluxes, and finally finds 
the minimum x 2 among all the redshifts for each model 
and among all the models. The output includes the best 
template, redshift, and scale factor to convert to absolute 
flux. These fits were only used to confirm the Hyperz re- 
sults because reddening was not taken into account in 
our code. 

4.2. FIR-Submillimeter 

The photometric measurements from the long- 
wavelength data were added to the previous optical and 
NIR photometry. The SED templates used previously 
cover the rest-frame wavelengths from 100 A to 3/im, so 
were further extended to the submillimeter with 4 long- 
wavelength templates (see Appendix A for a more de- 
tailed descript ion of the procedur e). The templates used 
were Arp 220 (Bressan et al. 2002) representing a ULIRG 
SED, M82 (Bress an et al.||2002p representing a star burst 
SEP, a nd 2 t emplates in the synthetic library from [Dale] 
et al. (2001) representing LIRGs (a=1.06) and quies- 
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TABLE 4 
OPTIC and NIR Photometric Data. 



# 


Name 


B 


&B 


V 




R 


an 


I 


CI 


z' 


°V 


J 




H 




K 




1 


F0023-0904 


24.44 


.16 


24.45 


.17 


23.41 


.15 


21.49 


.10 


21.69 


.11 


19.31 


.25 


18.83 


.16 






2 


F0129-0140 


23.69 


.15 


23.09 


.11 


22.68 


.04 


21.93 


.21 


21.99 


.08 


20.84 


.27 


19.23 


.15 






3 


F0152-0029 


>25.25 




>24.75 




>24.17 




>23.51 




>23.73 








>20.05 




19.34 


.12 


4 


F0152+0052 


>25.07 




>24.87 




>25.11 




>24.99 




>24.84 




>18.47 




>19.00 








5 


F0202-0021 


25.45 


.15 


22.61 


.07 


21.64 


.06 


19.85 


.07 


20.05 


.25 


18.31 


.20 


18.15 


.24 


16.72 


.23 


6 


F0216+0038 


>24.73 




22.91 


.11 


21.78 


.06 


20.20 


.07 


19.74 


.11 


18.84 


.20 


17.65 


.24 


16.37 


.23 


7 


F0916+1134 






22.65 


.11 


22.50 


.06 


21.84 


.07 


21.33 


.23 






19.69 


.15 






8 


F0919+1007 


>25.11 




23.61 


.11 


22.66 


.06 


20.48 


.07 


19.54 


.25 






17.85 


.15 


17.17 


.09 


9 


F0938+2326 






21.73 


.11 


20.88 


.06 


20.56 


.08 


20.05 


.02 






19.19 


.15 






10 


F0939-0128 






>24.91 




>25.15 








>24.54 








20.98 


.20 






11 


F0942+1520 


>24.56 




>24.31 




>25.08 




>24.26 












20.83 


.17 


18.13 


.15 


12 


F0943-0327 


>25.05 




>24.89 




>25.14 




23.57 


.21 










>20.07 








13 


F0950+1619 










23.74 


.10 


22.49 


.09 


>24.55 








>20.00 








14 


F0952+2405 














20.72 


.07 










17.26 


.23 






15 


F0955+2951 






24.13 


.13 


23.07 


.07 


21.91 


.08 


21.43 


.10 






20.31 


.14 






16 


F0955+0113 






>24.88 




>25.16 




>25.02 




>24.54 








>20.02 








17 


F0956-0533 






>24.91 




>25.15 




>24.97 












>19.06 








18 


F0958+2721 














18.96 


.07 










17.85 


.14 






19 


F1000-0636 














>25.00 












>20.02 








20 


F1008-0605 


21.71 


.16 


20.08 


.12 


18.94 


.20 


17.98 


.18 










16.32 


.14 






21a 


F1010+2527N a 


>23.78 




>23.15 




23.73 


.16 


22.23 


.16 


22.08 


.18 






20.49 


.17 


18.55 


.22 


21b 


F1010+2527S a 


>23.78 


.23 


>23.15 


.11 


22.48 


.06 


21.00 


.07 


21.02 


.08 






20.02 


.14 


17.61 


.25 


22 


F1010+2727 


>24.22 




>24.47 




22.44 


.25 


21.74 


.07 


22.27 


.17 






19.00 


.18 


18.08 


.15 


23 


F1014+1438 


>25.18 




>25.63 




24.46 


.11 


22.81 


.07 


22.46 


.08 






19.88 


.17 


18.89 


.13 


24 


F1016+1513 










>25.28 




>25.14 




>24.58 








>20.11 








25 


F1024-0031 


























17.49 


.33 






26 


F1027+0520 






18.17 


.11 


17.33 


.06 


16.67 


.07 


16.78 


.08 






15.14 


.14 






27 


F1039+2602 


>25.14 




24.61 


.34 


23.85 


.45 


22.33 


.28 


22.29 


.44 






20.24 


.15 


18.33 


.04 


28 


F1040+2323 










>25.20 




>25.09 




23.27 


.13 






20.66 


.17 






29 


F1116+0235 


21.06 


.15 


20.50 


.12 


19.52 


.06 


18.68 


.07 


18.83 


.08 






16.58 


.14 






30 


F1133+0312 










>25.20 




>25.07 




>24.55 








>20.12 








31 


F1140+1316 


>25.14 




>25.74 




>25.20 




>25.04 




>24.55 
















32 


F1147+2647 


>24.99 




>24.89 




25.24 


.10 


22.50 


.20 


22.82 


.12 






20.66 


.16 






33 


F1155+2620 


>25.00 




>24.91 




23.88 


.19 


22.30 


.08 


21.98 


.09 






19.80 


.17 






34 


F1158+1716 


























>19.06 








35 


F1202+0654 


























>19.02 








36 


F1211+3616 






>25.72 




>25.98 




>25.60 












>20.13 








37 


F1215+3242 


20.36 


.25 


18.86 


.21 


17.98 


.26 


17.23 


.27 


16.33 


.19 






15.83 


.14 






38a 


F1217-0529E a 






>24.88 




24.66 


.09 


22.34 


.11 


22.82 


.09 






18.90 


.26 






38b 


F1217-0529W a 






>24.88 




24.83 


.17 


22.98 


.09 


22.95 


.14 






20.73 


.17 






38 


F1217-0529S a 






>24.88 




25.47 


.19 


23.58 


.21 


>24.05 








>20.93 








39 


F1217+3810 


























15.52 


.14 






40 


F1218-0625 






24.13 


.20 


23.27 


.20 


21.91 


.28 


22.14 


.42 






19.48 


.18 






41 


F1218-0716 






23.82 


.14 


















21.54 


.20 






42 


F1234+2001 






>24.88 




>23.20 




19.85 


.28 


19.58 


.20 






18.63 


.16 






43 


F1237+1141 


>24.02 




23.99 


.07 


23.27 


.06 


22.47 


.07 


22.38 


.04 






19.64 


.29 


18.68 


.14 


44 


b 1315+4438 






23.41 


.30 


23.37 


.20 


22.66 


.66 


22.83 


.48 






>20.04 




21.32 


.15 


45 


F1329+1748 






20.78 


.11 


20.32 


.06 


19.68 


.07 


















46 


F1355+3607 










>23.42 




>23.87 




>23.68 








>21.07 




>20.86 




47 


F1430+3557 


>25.18 




>24.14 




>23.52 




>23.14 




>23.59 








>19.89 




18.86 


.13 


48 


F1435-0029 


>25.47 




>25.86 




>26.11 




>25.00 




22.91 


.12 






19.87 


.15 






49 


F1445+2702 














20.69 


.07 


















50 


F1447+1217 


>24.89 




>24.70 




24.00 


.11 


22.60 


.09 


22.16 


.11 






>20.25 




20.20 


.07 


51 


F1451+0556 


19.80 


.15 


19.59 


.11 


19.01 


.06 


18.81 


.07 


19.07 


.15 






17.80 


.14 


17.06 


.20 


52a 


F1458+4319NW b 


>25.03 




>24.91 




24.73 


.14 


22.30 


.09 


21.99 


.12 






21.45 


.18 






52b 


F1458+4319SE b 


>25.03 




>24.91 




24.41 


.09 


21.38 


.10 


21.77 


.09 






20.19 


.16 






52c 


F1458+4319E b 


>25.03 












21.17 


.08 


21.28 


.09 






19.29 


.33 






53 


F1505+4457 


>25.33 


.15 


23.98 


.07 


21.90 


.11 


21.18 


.07 


20.94 


.08 






18.46 


.16 


18.24 


.28 


54 


F1524+5122 










>23.87 




>23.17 




>23.55 








>20.31 




>20.96 




55 


F1644+2554 










>25.20 




>25.13 




>24.55 








>19.89 








56 


F2217-0837 


23.14 


.20 


20.30 


.13 


20.10 


.07 


19.03 


.07 


19.16 


.08 


17.77 


.24 


17.31 


.14 


16.59 


.03 


57a 


F2217-0138E a 


>24.98 




24.48 


.16 


22.42 


.08 


20.59 


.07 


20.88 


.09 


20.75 


.28 


18.64 


.15 


17.48 


.10 


57b 


F2217-0138W a 


>24.98 




>24.89 




24.60 


.15 


22.32 


.09 


22.01 


.14 


22.00 


.50 


20.39 


.16 


19.29 


.34 
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FlG. 9. — SED fits of 4 objects in the FIRST-BNGS sample. For each object, the top panel shows the best-fit SED and the actual broadband 
photometry data points with 1-sigma error bars in the rest frame of the galaxy. The bottom panel shows the distribution of best % 2 for the 11 
templates at each redshift. 
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FlG. 10. — Same as Figure 9 for 4 additional FIRST-BNGS sources. 



cent (a=2.5) SEDs. The best fit SEDs with each long- 
wavelength SED template is shown in Figures 18-21. 

For the four objects with FIR flux measurements, all 
seem to exclude an Arp-220 (ULIRG). Also the three at 
high redshift {z > 1) seem to exclude a LIRG-type SED. 
The only constrained model (from multiple Spitzer detec- 
tions) at z=2.90 quite definitely shows a M-82 (starburst) 
type SED. 



4.3. The IR Hubble Diagram 
In Figure 22 we plot the K-b and Hubble D i agram 



with literature data compiled by Willott et al. (2003). 
These data sets include more powerful radio galaxies 
and quasars with spectroscopic redshifts compared to 
the FIRST-BNGS sample with photometric redshifts (in- 
cluding the long wavelength SED solutions from section 
4.2 and H-band fluxes converted to K-band when neces- 
sary according to the best-fit SED). The literature mag- 
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FlG. 11. — Same as Figure 9 for 4 additional FIRST-BNGS sources. 



nitudes have been corrected to a common physical aper- 
ture (63.9 kpc) and also accounted for expected strong 
emission line flux according to a prescription by |Jarvis| 
et al.| (2001). Our K-band fluxes were corrected to the 
same physical aperture according to the GALFIT mod- 
els. The best-fit second-order polynomial of the litera- 
ture data traces closely a passively evolving stellar pop- 
ulation that forme d at z^lO from [Bruzual Sz Charlotl 
( |2003| ) as noted by |Willott et al.| fl2003| ). The agreement 



between the literature and the FIRST-BNGS sample is 
reasonable considering the accuracy of photometric red- 
shifts relative to the spectroscopic redshifts of the liter- 
ature. Also it is probable that our data is incomplete at 
faint flux levels (there are about 20 FIRST-BNGS sources 
not detected in H or K) since we were only surveying 
in H-band to a particular depth (H~21). The FIRST- 
BNGS sample seems to follow the passively evolving pop- 
ulation out to at least z = 5. 
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FlG. 12. — Same as Figure 9 for 4 additional FIRST-BNGS sources. 



5. SPECTROSCOPY 

Throughout the initial and follow-up phases of this 
project, we had the opportunity to observe these objects 
with various spectrographs on Mauna Kea to confirm 
the photometric redshifts. These observations were chal- 
lenging because of the faint ness of the science targets, 
particularly for slit alignment and guiding over long in- 
tegration times. However, the acquisition problems were 
largely solved because of the accurate offsets (from the 



Pueo-KIR AO observations) from the nearby bright AO 
guide star. This also helped with tracking using the AO 
star on the guider. Table 7 lists the observed targets, 
dates, and observing mode used. 

Several of the brig htest NIR sources in the sample were 
observed with SpeX (iRayner et al. 2003 ) on IRTF for sev- 



eral hours each in low-resolution prism mode to get the 
widest spectral coverage (0.8-2.5 fim) an d best sensitiv- 
ity. Reductions were done with Spextool (Cushing et al. 
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FlG. 13. — Same as Figure 9 for 4 additional FIRST-BNGS sources. 



2004; Vacca et al.||2003 ) including wavelength, flux, and 
telluric calibrations. 

Some of the optically brightest objects were ob- 
served w ith the Gemini- Nor th Multiobject Spectrograph 
(GMOS; |Hook et al.| ( [2003D ) in long-slit, nod and shuf- 
fle mode. Reductions were done using the GMOS IRAF 
reduction package with the nod and shuffle routines to 
carry out wavelength and flux calibrations. No telluric 
correction was done. 



A few of the most interesting sources were ob- 
served with the OH-suppress i ng IR Imaging Spectro- 
graph (OSIRIS; |Larkin et ah] fl2006| )) on Keck II. This 
AO-fed integral field spectrograph utilizes a lenslet ar- 
ray to obtain full near-IR broadband spectral coverage 
(z, J, H or K) at R~4000 resolution at approximately 
1000 well-sampled spatial locations. It is more efficient 
than a traditional longslit spectrograph as no light from 
extended objects is blocked out by the slit. In addi- 
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FlG. 14. — Same as Figure 9 for 4 additional FIRST-BNGS sources. 



tion, the high resolution enables observations between 
the bright OH sky lines, resulting in lower background 
noise. The data reductions were done with the OSIRIS 
data reduction pipeline written in IDL, which provided a 
wavelength rectified data cube for each object observed. 
Flux calibration, telluric correction, and spectral extrac- 
tion were accomplished with our own IDL scripts. 

For each object observed with each instrument, we fit 
the flux and wavelength calibrated 1-D spectra to the 



same SED library of 7 templates we used for the pho- 
tometric redshifts, allowing redshift and luminosity to 
vary, and matching the spectral resolution to the data 
via convolution. We used IDL-based routines to perform 
a least-squares fit of the templates to each spectrum. 
Figure 23 shows examples of the best fits of template 
models with the data and Table 8 lists the results from 
these fits. A number of real and spurious features were 
found by examining the individual beam exposures. We 
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FlG. 15. — Same as Figure 9 for 4 additional FIRST-BNGS sources. 

discuss each object briefly below. 

F0023-0904 with GMOS has a confirmed [Oil] 3727 
detection at 7280A and forbidden [Nelll] 3869 at 7550A. 

F0916+1134 with GMOS has a confirmed [Oil] 3727 
detection at 6640A , and contamination from a nearby 
source is present at 8160A. 

F1014+1438 with GMOS shows a steep break at about 
6 100 A which was well fit to the Lyman break of a SB2 
SED.This is consistent with the Hyperz fit. 



Fl 234 +2001 with GMOS has a steep break at 7900A. 
This was successfully fit to the Lyman break, but with 
no other features is not conclusive. The two OSIRIS 
spectra show faint continuum with no large features, and 
the SpeX data shows a faint continuum and possibly the 
Lyman break just at the edge of sensitivity. All four fits 
roughly agree are also consistent with the photometric 
Hyperz fit of a z~5.5 starburst. 

F 1237+ 1141 with GMOS is featureless and flat ex- 
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FlG. 16. — Same as Figure 9 for 4 additional FIRST-BNGS sources. 

cept some bad sky subtraction at 7900A to 8000A. The 
OSIRIS data shows a break at 1.6 jam and perhaps [Oil] 
3727 at 1.5/im. 

F U5 1+0556 with GMOS shows the edge of [CIV] 1550 
at 5400A and [CIII] 1909 at 6800A. SpeX shows [Mgll] 
2798 at 1.0/im, [Oil] 3727 at 1.35/xm, H/3/[0III] 5007 at 
1.73/xm, and Ha at 2.32/im. The bad fit to H/3 is proba- 
bly because in the templates there is a sharp transition 
in Kf3 flux between SB2 and SB3 templates. OSIRIS 



also shows H/3 at 1.73/im, [OIII] at 1.79/im, and Ha at 
2.32/im. The equivalent widths of Ha and H/3 were mea- 
sured by Gaussian fitting and given in Table 9. If we 
adopt the best S/N (SpeX) measurements, the derived 
reddening is E(B-V)=1.73. 

Fl 505+4457 with GMOS shows the Balmer break fit 
at about 8000A. The SpeX data is consistent with the 
Balmer break and a smooth continuum. 

F2217-0138 with GMOS shows a break which was fit 
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FlG. 17. — Same as Figure 9 for 3 additional FIRST-BNGS sources. 

to the Lyman break. 

The rest of the spectroscopy did not produce any well- 
constrained redshifts. A few objects were detected in 
continuum at relatively low S/N, but without a definitive 
spectral break or emission lines, the fits could only con- 
firm that the photometric redshifts were consistent with 
the spectra. We present these less-constrained spectro- 
scopic results in Table 10. 

The well-constrained spectroscopic redshifts may be 



compared to the photometric redshifts derived in the pre- 
vious section. Figure 24 shows the comparison between 
the two. 

The photometric redshifts match the spectroscopic 
ones well with a RMS of <5z/(l+z)=0.146. This con- 
firms that the photometric redshifts, though less accu- 
rate than spectroscopic ones, are still valid overall. His- 
torically, Monte Carlo simulations for photometric red- 
shifts from optical+NIR surveys have determined a typi- 
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Best-fit SEDs for Fl 116+0235 with ULIRG (solid), LIRG (dotted), starburst (dot-dashed), and quiescent (dashed) long wavelength 



cal error of £z/(l+z)=0.08 (Labbe et al.|2003| ), but more 
conservatively estimate 0.1. This is an encouraging re- 
sult considering our spectral template library may not 
be an accurate representation of the high redshift stellar 
populations we observed. It also adds confidence to our 
photometry measurements that were made in difficult re- 
gions (near a bright star with diffraction spikes and near 
the edge of the detector). 

6. VLA OBSERVATIONS 

VLA observations were carried out on Nov 28th, 2004, 
Feb 28th, 2006, and April 20th, 2006 with the VLA in A- 
array configuration at 3.6cm. In addition, a VLA archive 
search resulted in 2 observations of F1445+2702 (known 
as IRAS14434+2714) on April 13th 1998 and January 
14th, 2001 with the same VLA configuration. We used 
calibrators 0016-002, 0137+331, 0219+013, 0954+177, 
0956+252, 1331+305, 1335+457, 1415+133, 1436+233, 
1500+478, 2229-085 for flux and phase calibration. All 
data were reduced using standard packages within the 
Astronomical Image Processing System (AIPS). Some 
sources were observed multiple times and were coadded 
to produce the final map. 

With the VLA A-array at 3.6cm, the largest angular 
scale detectable is 7.0 arcseconds and the primary beam 
is about 0.24 arcseconds. Flux densities were measured 
using the IRAF "phot" package and converted to flux 
units. Background noise was also measured to derive 
flux errors and upper-limits for non-detection. For mul- 
tiple sources (double, triple, distorted morphologies), the 
fluxes were added together. Figure 25 shows example ra- 
dio maps from FIRST and from our VLA A-array obser- 
vations. Table 11 shows our measurements along with 



the best-fit photometric or spectroscopic redshifts plus 
spectral indices. The spectral index, a, where S v oc 
was calculated for each discrete radio source in the field 
as well as the total flux of all the sources coadded. 

Plotting spectral index of each object's total fluxes, a 
versus 20-cm flux, redshift, absolute magnitude, there is 
no clear indication of any trends with a in this data set 
(Figure 26) though there are more steep-spectrum (a < 
—0.5) sources observed than flat-spectrum (a > —0.5) 
sources (9 vs. 3). 

7. DISCUSSION OF INDIVIDUAL OBJECTS 

Below is a discussion the properties of the 58 individual 
sources in the FIRST-BNGS sample. 

F0023-0904 was detected in all bands, giving a pho- 
tometric redshift of 1.49 with an Sc SED. Spectroscopy 
revised the redshift to 0.946 with an [Oil] emission detec- 
tion. The FIRST radio peak is offset about 3" from the 
galaxy, but the outer contours are elongated in the di- 
rection of the galaxy, suggesting a jet structure. VLA X- 
band morphology was inconclusive, but provided a spec- 
tral slope of -1.04. 

F0129-0140 had 2 solutions for photometric redshifts 
at 0.65 with an SB2 SED and 2.44 with an SB4 SED. 
Without any discernible break, it was not obvious which 
is the better solution. Additionally, spectroscopy gives a 
rough redshift of 1.7 from the spectral slope in the 5500- 
8000 Angstrom range which does not break the degen- 
eracy. A compact moderate radio flux with faint optical 
flux suggests a higher redshift as the low redshift would 
give such a low luminosity (9* 10 9 L snn ), it would be dif- 
ficult to reconcile with an AGN host. We will use the 
high photometric redshift solution because it was more 
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constrained than the spectroscopic fit. 

F0152-0029 was undetected in all bands except K. It 
does not have a constrained redshift so is either at z>8.0 
where the Lyman break is redshifted out of z-band (see 
figure 20) or a highly reddened object. The radio source 
is compact with a relatively high radio 20cm luminosity 
(25mJy) and steep spectral index (-1.09), suggesting an 
AGN source. 

F0152+0052 was not detected in any band (B, V, R, 
I, z', J, and H). A triple morphology in X-band, and 
high 20cm flux (18mJy) and steep spectral index (-0.73) 
suggests an high redshift AGN. 

F 020 2- 00 21 appears to be well fit by a low redshift el- 
liptical SED (z = 0.58), but also has a secondary solution 
at z = 4.2. The brightness and redshift inferred from the 
K-z relation derived in section 2.4 (z ~ 0.9) supported 
the first solution and was roughly confirmed with spec- 
troscopy (z = 0.6). An elongated radio morphology and 
lower radio 20cm flux (3mJy) with a steep spectrum also 
is consistent with a giant elliptical with an AGN at low 
redshift. 

F 021 6+ 00 38 has a SED consistent with the SO tem- 
plate at z = 0.65 which is lower than that inferred from 
the K-z relation derived in section 2.4 (z ~ 0.9), but 
was roughly confirmed with spectroscopy. A moderate 
20cm flux (12mJy) and double radio lobe morphology 
with steep spectral index (-1.39) suggests a relatively 
powerful AGN, perhaps near the end of an active ac- 



cretion stage. 

F0916+1134 appears to be a barred-spiral from the 
residuals in surface brightness fitting (Stalder & Cham- 
bers, in prep.). K-band photometry was excluded from 
photo-z fitting (due to the poor residuals in the profile 
fit), which found a best redshift of 1.12 with a starburst 
SED. A confirmed [OIII] 3727A detection with spec- 
troscopy gives a final redshift of 0.78. 

F0919+1007 appears to have a low redshift (z = 0.79) 
SO SED. Moderate radio 20cm flux (9m Jy) with elon- 
gated morphology is also consistent with a low-redshift 
giant elliptical with an AGN. 

F0938+2326 has 2 solutions for photometric redshifts 
at 0.8 with a SB2 SED and 3.88 with SB3 SED. The mod- 
erate to high 20-cm flux (8.1mJy) with compact morphol- 
ogy, along with the relatively faint optical flux, favors the 
higher redshift solution. 

F0939-0128 was not detected in V, R, z', or H bands. 

F0942+1520 was undetected in all bands except H and 
K. Because of the large H-K color, Hyperz could fit the 
break with a redshift of 3 to 6.5. The high radio 20cm 
flux (12mJy) and double- lobe morphology is consistent 
with a high-redshift AGN. 

F 09 43-0321 was undetected in B, V, R, and H, but 
barely detected in I (< 5a) and has a very high radio 
20cm flux (99mJy) and double-lobe morphology. 

F 0950 +1619 was detected in R and I though the ID is 
slightly offset (about 5") from the radio peak so it is pos- 
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(dotted), starburst (dot-dashed), and quiescent (dashed) long wavelength templates. 



sibly another blank field. There were also no detections 
in z and H bands. 

F0952+2405 was detected in both I and H bands. It 
is bright in H band, and a faint radio source (1.3mJy). 
It is probably not at high redshift (the K-z relation gives 
likely redshift of z=0.76). 

F0955+2951 has a distorted radio morphology with 
the radio peak offset from the optical ID by about 5". 
The optical and IR morphologies appeared extended. 
The best fit photometric redshift is a z=4.4 starburst. 
Since it is not likely the true optical ID for this radio 
source, we will not use it in the subsequent analysis. 

F0955+0113 was not detected in any observed filter 
(V, R, I, z' and H). 

F0956-0533 was not detected in V, R, I, or H bands. 

F 095 8 +27 21 similar to F0952+2405, detected in I and 
H and relatively bright with faint 20cm radio source 
(2.4mJy), so probably not at high redshift. 

F1000-0636 was not detected in I or H filters. 

F1008-0605 was detected and relatively bright in all 
observed filters (B, V, R, I, and H). A break was found 
and fit by Hyperz to either a z = 0.27 elliptical SED 
or z = 3.9 SB4 SED with little extinction. If the high 
redshift solution were correct, it would be hard to rec- 
oncile this particular object with a hierarchical merging 
scenario especially since it would have an extremely lumi- 
nous stellar population (> 1O 13 L ) dominating its rest 



frame UV light, which has to be assembled in less than 
2.3Gyr (age of the universe at that redshift). We there- 
fore adopt the low redshift solution. 

F1010+2527 has a close double morphology in most 
bands. There was not any observed evidence of interac- 
tion, but the photometric redshifts of the two galaxies 
were similar (z = 4.41,z = 4.63), which suggests they 
are associated. Both are consistent with starburst SEDs. 
The south source's parameters were used in the subse- 
quent analyses because the fit was better constrained. 
The unresolved, low 20cm flux (lmJy) radio source also 
supports that this is not an aligned object. 

Fl 01 0+2727 has an SED consistent with a high red- 
shift Sb model (z = 4.53). A moderate 20cm radio source 
(6mJy total) and double-lobe morphology also suggests 
a high redshift. 

F1014+1438 appears to be a starburst at high redshift 
(z = 3.95), but a low redshift (z = 0.5) solution was also 
found. Because of the faintness of the source and the 
presence of a steep break observed with spectroscopy, we 
chose the high redshift solution. Double-lobe moderate 
20cm flux (33mJy integrated) radio source is also consis- 
tent with a high redshift AGN. 

F1016+1513 was undetected in R, I, z', and H bands. 

F1024-0031 was detected in H, and with the most pow- 
erful radio 20cm flux (158mJy) in the sample, is almost 
certainly an AGN. Based on the H-band brightness, it is 
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FlG. 21. — Best-fit SEDs for F2217-0138E with ULIRG (solid), LIRG (dotted), starburst (dot-dashed), and quiescent (dashed) long wavelength 
templates. 




z 

FlG. 22. — The K-z relation for 3CRR, 6CE, 6C*, and 7CRS radio galaxies from [Willott et al.N2003) (stars) with spectroscopic redshifts and 
the FIRST-BNGS sa mple with photometric redshifts (boxes), all corrected to a physical aperture of 63.9 kpc. A pol ynomial fit f solid line) to t he 
| Willott et al.| |2003) data is K = 17.34 + 4.385 * log(z) - 0.4439 * log(z) 2 , and single stellar population models from |Bruzual fe Charlotj \200^ is 
shown - The models are for a non-evolving population (dot-dashed) and passively evolving populations from Zf = 10 (short dashed) and ~z~f~= 5 
(long dashed). 
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FlG. 23. — Examples of spectroscopic fitting of 5 FIRST-BNGS objects taken with various instruments (solid) with best fit template model 
(dashed). The object name, template model, and redshift are indicated in each panel. First row: F0023-0904 with GMOS (left) and F0916+0038 
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FlG. 24. — Comparison of redshifts derived from spectroscopic fitting and photometric fitting. Bottom panel is the residuals between spectroscopic 
and photometric redshifts. Calculated RMS <5z/(l+z)=0.146. 



probably at a redshift around 0.8. 

F1027+0520 was fit to a z=0.6 starburst SED tem- 
plate. It is a bright galaxy, so is potentially at lower 
redshift due to our low confidence in our photometric 
measurements for bright extended objects and determin- 
ing a good sky level around them. It has a compact radio 
source with relatively high 20cm flux (23mJy). 

F 1039+2602 has an SED best fit by a z = 3.62 young 
starburst SB3. Its moderate radio 20cm flux (12mJy) 
and double morphology is also consistent with high red- 
shift. 

F1040+2323 has a weak (1.6mJy) radio source blank 
in R and I, but faintly detected in z' and H. 

F1116+0235 has a best-fit photometric redshift at 
0.62 which is consistent with a continuum fit from 
spectroscopy. It also has a 2a detection at 850/im 
from SCUBA, which can only exclude a Arp 220-type 



(ULIRG) SED at the same redshift as the optical/NIR 
photometric redshift. 

Fl 133+0312 was undetected at R, I, z, and H bands 
and is at the lower threshold of the FIRST survey 
(0.8m Jy). 

Fl 140+1316 was not detected at B, V, R, I, and z' A 
double- lobed morphology suggests an AGN. 

F1147+2647 has a possible faint optical ID about 6 
arcseconds to the west of the weak radio peak (0.7mJy), 
so it is probably not the proper ID and will not be used 
in later analyses, though the best fit SED is a z=4.75 
SB1 type SED. 

F1155+2620 has a lmJy radio source with an opti- 
cal counterpart about 3 arcseconds to the southeast of 
the radio peak. Hyperz gives the best-fit photometric 
redshift at 4.5 with a starburst SED. 

F1158+1716 undetected in H. 
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FlG. 26. — Distribution of en with redshift (left), radio flux (center), and V-band absolute magnitude (right). Shapes represent best-fit optical/NIR 
SEDs from Hyperz, triangles are early-type (E and SO), boxes are late-type (Sa,Sb,Sc), and pentagons are starbursts SB1-SB6. The horizontal line 
show the division between steep (a <-0.5) and flat (a >-0.5) spectrum types. 
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TABLE 5 
Hyperz Results 



Name 


Zphot 






M 




SED 


F0023-0904 


1 


49 


+o 

-0 


14 
07 


-24 


10 


Sc 


F0129-0140 


2 


44 


+o 

-0 


40 
90 


-24 


73 


SB4 


F0202-0021 





58 


+o 

-0 


08 
02 


-22 


13 


E 


F0216+0038 





65 


+o 

-0 


30 
11 


-22 


09 


SO 


F0916+1134 


1 


12 


+o 

-0 


07 

05 


-22 


00 


SB3 


F0919+1007 





78 


+o 

-0 


06 
03 


-23 


08 


SO 


F0938+2326 


3 


88 


+o 

-0 


14 
07 


-26 


42 


SB3 


F0942+1520 a 


3 


36 


+ 3 
-0 


50 
54 


-26 


45 


SO 


F0955+2951 


4 


39 


+o 

-0 


14 
13 


-25 


77 


SB4 


F1008-0605 


n 
u 




+o 

-0 


15 

08 


-21 


32 


E 


F1010+2527N 




41 


+o 

-0 


33 
38 


-26 


57 


SB4 


F1010+2527S 


A 


Do 


+o 

-0 


08 
06 


-27 


32 


SB2 


F1010+2727 




53 


+o 

-0 


20 
11 


-27 


87 


Sb 


F1014+1438 


4 


47 


+o 

-0 


34 
22 


-26 


58 


SB5 


F1027+0520 





60 


+o 

-0 


07 
04 


-25 


18 


SB1 


F1039+2602 


3 


62 


+o 

-0 


69 
26 


-26 


18 


SB3 


F1116+0235 





62 


+o 

-0 


04 
04 


-23 


20 


SB2 


F1147+2647 


4 


75 


+o 

-0 


15 
11 


-25 


98 


SB1 


F1155+2620 


4 


54 


+o 

-0 


69 
49 


-26 


79 


SB3 


F1215+3242 





1 1 


+o 

-0 


29 
11 


-19 


86 


SO 


F1217-0529E 


4 


95 


+o 

-0 


21 

06 


-29 


93 


so 


F1217-0529W 




97 


+o 

-0 


32 
20 


-28 


96 


E 


F1217-0529S 


4 


82 


+o 

-0 


71 

80 


-23 


50 


SB2 


F1218-0625 


3 


98 


+o 

-3 


53 
99 


-27 


25 


Sb 


F1234+2001 


5 


40 


+o 

-0 


70 
41 


-27 


35 


SB4 


F1237+1141 


2 


38 


+o 

— 


19 

14 


-24 


34 


SB2 


F1315+4438 


2 


77 


+o 

-0 


16 

28 


-22 


49 


SB2 


F1447+1217 


4 


70 


+o 

-0 


13 
18 


-25 


11 


SB3 


F1451+0556 


2 


71 


+o 

-0 


15 
05 


-26 


70 


SB2 


F1458+4319NW 


5 


05 


+o 

-0 


25 
10 


-24 


40 


SB2 


F1458+4319SE 


5 


03 


+o 

-0 


34 
19 


-26 


06 


SB2 


F1505+4457 





55 


+o 

-0 


10 

08 


-20 


52 


Sa 


F2217-0837 





26 


+o 

-0 


04 
02 


-20 


16 


E 


F2217-0138E 


4 


55 


+o 

-0 


06 
21 


-27 


81 


SB2 


F2217-0138W 


4 


99 


+o 

-0 


in 

25 


-26 


35 


SB2 



a Only 2 bands detected, Lyman break between H and K. 

F1202+0654 undetected in H. 

F 121 1+3616 undetected in V, R, I, and H. 

F 1215+4342 appears to be a bright SO type galaxy 
at z=0.1 with a relatively high 20cm flux (45mJy) radio 
source. Although there was a solution at high redshift, 
the extreme optical luminosity (> 10 13 L snn ) estimated 
for that redshift and extended radio morphology suggest 
the low redshift solution is more likely. 

F 1211-0529 is another optical double (possibly triple) 
source which also has similar photometric redshifts 
(z=4.95 and z=4.97 for the east and west sources respec- 
tively). However, the radio peak is offset slightly (about 
3") to the south, closer to the faint (R=25.5) south op- 
tical source. 

F1217+3810 was detected at H-band and is relatively 
bright, so it is probably at low redshift (z ~ 0.3) with a 
weak (0.8mJy) radio source. 

F1218-0625 was detected in V, R, I, z' and H. Two 
photometric redshift solutions were found though any 
redshift between and 4.5 gives reasonable fits. The 
compact moderate 20-cm flux (4.3mJy) and faint optical 
flux suggest the higher redshift is more likely, but with a 
large possible range. 

F1218-0716 was barely detected in V and H bands. It 
also is a weak (0.8mJy) and distorted radio source. 



F 1234 +2001 was identified about 4" to the southwest 
from the moderate radio peak (5mJy) of the FIRST 
source J123432. 9+200134. This seems a bit far given the 
position accuracy of FIRST (about 1") and our imaging 
data (about 0.3") derived from the USNO-A2.0 catalog. 
Since the only optical/IR source in the field around the 
radio source is F1234+2001 there are 3 possibilities: 1) 
The source we have identified is unrelated or a compan- 
ion to the radio source J123432. 9+200134 in which case 
the radio source host is below the detection threshold 
of H=19.74 (Sa); 2) The optical emission is offset from 
the c enter of the host gal axy due to the alignment ef- 
fect ([Chambers et al.||1987[) though the radio source is 



too weak to be regarded as a powerful radio source at 
any epoch which makes this scenario unlikely; or 3) The 
radio source is intrinsically asymmetric due to relativis- 
tic beaming, so the radio centroid of J123432. 9+200134 
is not centered on the host galaxy F 1234+2001. For the 
remainder of the paper we accept that optical/IR source, 
F1234+2001, to be the host of J123432.9+200134. The 
confidence in this ID is strengthened from the subse- 
quently derived redshifts. 

It should be noted that a diffraction spike passes 
through the optical ID in V, R, I, and z' . GALFIT 
successfully modeled this and the good residuals raise 
our confidence in the photometric measurements. The 
best fit SED to the broadband photometry is a z = 5.4 
starburst mainly from fitting the Lyman break between 
R and I bands. Though the R-band image is relatively 
shallow, the V-band imaging is deeper which did con- 
firm that this object was at least z > 3. At z = 5.40, 
the imaging data spans the UV wavelength range in the 
galaxy's rest-frame from 800A (V-band) to 2800A (H- 
band). The galaxy is unresolved in all bands, including 
the AO H-band (0.29" FWHM) making its physical ex- 
tent on order or less than 1.7 kpc (1.0" corresponds to 
6.0 kpc at z ~ 5.5). 

Unfortunately if the observed spectral break is the Ly- 
man break, there are few powerful emission lines to ob- 
serve spectrally in this wavelength range. The Lyman 
break would be observed at about 8000A which will make 
it the primary spectroscopic feature for our fitting rou- 
tine. Because it is so bright (1=19.85), absorption fea- 
tures may be also detectable with a 8-10 meter class red- 
sensitive spectrograph. 

The GMOS spectroscopy does show a strong break 
around 7900A and confirms our photometric redshift fit 
that the Lyman break is between R and I bands. The 
signal to noise ratio of the spectrum was insufficient to 
identify any absorption features. 

The SpeX spectroscopy also shows the Lyman break 
though just at the edge of sensitivity. Both spectra are 
also consistent with the photometric data, give similar 
redshifts and are within the error of the photometry- 
derived redshift (Tables 5 and 8). A weighted aver- 
age of the three redshifts gives a best estimate at z = 
5.53 + 0.06. These three independent measures suggest 
that if the optical ID, F1234+2001, is associated with 
the radio source, J123432. 9+200134, it would make it 
the most distant known radio galaxy. However, this red- 
shift should be confirmed with deeper spectroscopy. 

F1237+1141 is consistent with a moderately high red- 
shift (z ~ 2.5) starburst SED. The radio morphology is 
a triple system at low radio 20cm flux (lmJy) with the 
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TABLE 6 

Spectroscopy Results 



Name 


Instruments 


%spec 




M v 


SED 


F0023-0904 a 


GMOS 


0.95 


0.01 


-19.41 


SB2 


F0916+1134 a 


GMOS 


0.78 


0.01 


-20.53 


SB5 


F1014+1438 b 


GMOS 


3.95 


0.15 


-23.80 


SB2 


F1234+2001 b 


GMOS,SpeX,OSIRIS 


5.56 


0.07 


-27.59 


SB4 


F1237+1141 b 


GMOS, OSIRIS 


3.01 


0.04 


-23.42 


SB3 


F1451+0556 c 


GMOS,SpeX,OSIRIS 


2.57 


0.01 


-26.76 


SB5 


F1505+4457 b 


GMOS,SpeX 


1.02 


0.03 


-20.45 


Sa 


F2217-0138E b 


GMOS 


4.66 


0.26 


-27.99 


SB2 



[OI/]3727 detected 13 break detected [MgII]2799, [OJ/]3727, Hp, [O///]5009, and H a detected (possibly [C///1909) 



TABLE 7 
F1451+0556 Equivalent Widths 









H(3 




Instrument 


EW 




EW 


O- EW 


SpeX 


1515A 


519A 


74. 9A 


53. 6A 


OSIRIS 


2266A 


3374A 


150A 


36A 



ID corresponding with the center source consistent with 
this redshift. It was also observed with both IRAC and 
MIPS, which changed the best fit redshift to 2.9 to re- 
flect the 5.8/im data point, but with large error bars. 
The SED beautifully fits a z = 2.9 starburst with a long 
wavelength starburst SED (M82) from a rest wavelength 
of 110 nm to 40 microns. This is a superb illustration 
of the potential of Spitzer for both photometric redshifts 
and studies of the stellar environment of these objects. 
Spectroscopy better constrained the redshift to 3.01. 

F1315+4438 is the faintest detected object in the 
HzRG candidate sample (K = 21.32). Hyperz found 
2 photometric redshifts, a SB1 SED at z = 1.1 and a 
SB2 SED at z = 2.8. The extremely faint optical flux 
favors the high redshift solution though a low radio 20cm 
flux (0.9mJy) with a slightly distorted morphology may 
indicate that the radio emission is from star formation. 
However, this interpretation would support an even lower 
redshift than z ~ 1. 

F 1329 +1748 is bright in V, R, and I with a weak 
(1.3mJy) radio source. 

F1355+3607 is undetected with 5a limit at K > 20.86. 
A distorted, faint (3mJy) radio source is consistent with 
a low redshift starburst, perhaps highly reddened. 

F1430+3557 was undetected in all bands except K so 
it does not have a photometric redshift until adding the 
MIPS 24/im and 160/im Spitzer archive data, the best-fit 
SED is a starburst at z = 3.5 with starburst or quiescent 
long wavelength SED (not LIRG or ULIRG). The radio 
source 20cm flux is moderate (9mJy), together with the 
lack of resolution of the radio source is also consistent 
with an AGN at high redshift. 

F1435-0029 was undetected in B, V, R, and I, but 
detected in z' and H. It has a moderate radio source 
(lOmJy), and an apparent significant break between I 
and z', so probably at high redshift (z ~ 3.6). 

F 1447+ 1217 is probably at very high redshift (z = 
4.70). It is one of the most powerful radio source (49mJy) 
in the sample, but a good SED fit, suggests little QSO 
contamination, so we might consider objects at lower ra- 
dio 20cm flux also safe from contamination. 



F1451+0556 is unresolved and has a peculiar color 
(R — K = 0.27), which most resembles a young star- 
burst SB2 SED at z=2.71, though not a very good fit. 
Spectroscopy confirms a redshift of 2.567 with several 
emission lines. This suggests significant AGN influence 
on the SED of this galaxy explaining the peculiar color 
and compactness. 

F1458+4319 has a distorted radio source making it 
difficult to identify an optical counterpart. The brighter 
two sources have photometric redshifts at z=5.1 and 5.0 
with starburst SEDs. The extreme I-R break seems con- 
sistent with the SED fit. 

F1505+4457 has an SED consistent with an Sa galaxy 
at low redshift (z = 0.55) though spectroscopy data sug- 
gests a higher redshift (z=1.02). The optical ID lies be- 
tween the 2 sources (which are at different fluxes. This 
seems to suggest a radio jet with the axis almost per- 
pendicular to the plane of the sky with the north lobe 
pointed away. With a total 20-cm radio of about 17mJy, 
with double-lobe morphology, it is consistent with this 
moderate redshift. 

F1524 +5122 is undetected with 5a limit at K > 20.96. 
The relatively high radio luminosity (22mJy) suggests 
an AGN. It is potentially a very high redshift object, 
possibly the host galaxy could be highly reddened. 

F1644+2554 was undetected at R, I, z', and H also 
was not detected with SCUBA so is not constrained as 
it has not been detected in any band. 

F2217-0837 has 2 Hyperz solutions, a z = 0.3 ellip- 
tical and z = 4.1 SB2 SED. The low redshift solution 
seems more likely as it is would otherwise be an ex- 
tremely bright (> 10 13 L snn ) galaxy at high redshift. A 
weak (lmJy), elongated radio source also supports the 
low redshift solution. 

F2217-0138 has two possible optical IDs, a brighter 
source about 1" to the east of the radio peak with a 
photometric redshift of z = 4.6 (confirmed with spec- 
troscopy), and a fainter source 1" to the west with a pho- 
tometric redshift z = 5.0. An upper limit from SHARC2 
does not help the optical ID nor constrain any of SED 
fitting. A powerful (26mJy) radio source suggests an 
AGN. 

F2354-0055 was undetected in B, V, R, I, z' and H 
filters, and a compact, faint, steep spectrum radio source 
is probably from an AGN, perhaps at high redshift. 

8. SUMMARY 

A set of 58 VLA FIRST survey sources that lie within 
the isoplanatic patch of a bright natural guide star 
(BNGS) was constructed to search for high redshift ra- 
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TABLE 8 
Other Spectroscopy 



Name 


Instruments 


%spec 




M v 


SED 


F0129-0140 


GMOS 


1.71 


0.40 


-22.80 


SB5 


F0202-0021 


GMOS,SpeX 


0.60 


0.30 


-22.15 


E 


F0216+0038 


GMOS,SpeX,OSIRIS 


0.56 


0.30 


-21.49 


SO 


F0919+1007 


OSIRIS 


0.79 


0.30 


-23.08 


SB5 


F1116+0235 


SpeX 


0.31 


0.20 


-20.65 


Sc 


F1447+1217 


GMOS 


4.60 


0.50 


-24.79 


SB5 


F2217-0837 


SpeX 


0.31 


0.20 


-20.16 


E 



TABLE 9 

VLA 3.6cm and 20cm Fluxes and 3-cr upper-limits. 



Name 


z 




M v 


SED 


S8.3 


C"8.3 


Si. 4 


or. 4 


a 


morph a 


F0023-0904 


0.946 


+ 0.005 
-0.005 


-19.41 


SB2 


2.56 


0.09 


15.28 


0.13 


-1.04 


elong 


F0129-0140 


2.76 


+ 0.27 
-0.44 


-25.07 


SB6 


0.83 


0.08 


3.70 


0.14 


-0.87 




F0152-0029 








3.79 


0.09 


24.69 


0.17 


-1.09 


unres 


F0152+0052E 










2.17 


0.09 


9.88 


0.17 


-0.88 


tpl 


F0152+0052W 










1.92 


0.09 


6.97 


0.17 


-0.75 


tpl 


F0152+0052C 










0.73 


0.09 


<0.6 




>0.11 


tpl 


F0152+0052tot 










4.83 


0.09 


16.85 


0.17 


-0.73 


tpl 


F0202-0021 


0.58 


-0.02 


-22.13 


E 


<0.27 




2.79 


0.11 


<-1.36 


elong 


F0216+0038E 


0.65 


+ 0.30 
-0.11 


-22.09 


SO 


2.32 


0.08 


17.84 


0.10 


-1.19 


dbl 


F0216+0038W 


0.65 


+ 0.30 
-0.11 


-22.09 


SO 


0.42 


0.08 


12.10 


0.10 


-1.96 


dbl 


F0216+0038tot 


0.65 


+ 0.30 
— 11 


-22.09 


so 


2.74 


0.08 


29.94 


0.10 


-1.39 


dbl 


F0916+1134E 


0.78 


+ 0.01 
-0.01 


-20.53 


SB5 


1.38 


0.20 


3.45 


0.15 


-0.53 


elong 


F0916+1134W 


0.78 


+ 0.01 
-0.01 


-20.53 


SB5 


<0.6 




1.09 


0.15 


<-0.35 


elong 


F0916+1134tot 


0.78 


+ 0.01 
-0.01 


-20.53 


SB5 


1.38 


0.20 


4.54 


0.15 


-0.69 


elong 


F0919+1007 


0.79 


+ 0.03 
-0.03 


-23.08 


SB5 


9.71 


0.25 


8.97 


0.15 


0.05 


elong 


F0942+1520E 


3.36 


+ 3.56 
-0.54 






6.17 


0.16 


11.70 


0.17 


-0.37 


dbl 


F0942+1520W 


3.36 


+ 3.56 
-0.54 






0.72 


0.16 


4.00 


0.17 


-1.00 


dbl 


F0942+1520tot 


3.36 


+ 3.56 
-0.54 






6.89 


0.16 


15.70 


0.17 


-0.48 


dbl 


F1010+2527 


4.63 


+ 0.08 
-0.06 


-27.32 


SB2 


<0.33 




1.31 


0.15 


<-0.80 


unres 


F1010+2727N 


4.53 


+ 0.20 
-0.11 


-27.87 


Sb 


<0.33 




2.36 


0.15 


<-1.15 


dbl 


F1010+2727S 


4.53 


+ 0.20 
-0.11 


-27.87 


Sb 


<0.33 




3.58 


0.15 


<-1.39 


dbl 


F1010+2727tot 


4.53 


+ 0.20 
-0.11 


-27.87 


Sb 


<0.33 




5.94 


0.15 


<-1.69 


dbl 


F1014+1438E 


3.95 


+ 0.15 
-0.15 


-23.80 


SB2 


<0.33 




17.14 


0.17 


<-2.30 


tpl 


F1014+1438W 


3.95 


+ 0.15 
-0.15 


-23.80 


SB2 


<0.33 




15.60 


0.17 


<-2.25 


tpl 


F1014+1438C 


3.95 


+ 0.15 
-0.15 


-23.80 


SB2 


4.11 


0.11 


<0.6 




>1.12 


tpl 


F1014+1438tot 


3.95 


+ 0.15 
-0.15 


-23.80 


SB2 


4.11 


0.11 


32.74 


0.17 


-1.21 


tpl 


F1237+1141N 


3.01 


+ 0.04 
-0.04 


-23.42 


SB3 


3.92 


0.05 


11.91 


0.17 


-0.65 


tpl 


F1237+1141S 


3.01 


+ 0.04 
-0.04 


-23.42 


SB3 


0.45 


0.05 


4.56 


0.17 


-1.35 


tpl 


F1237+1141C 


3.01 


+ 0.04 
-0.04 


-23.42 


SB3 


1.06 


0.05 


1.36 


0.17 


-0.15 


tpl 


F1237+1141tot 


3.01 


+ 0.04 
-0.04 


-23.42 


SB3 


5.43 


0.05 


17.83 


0.17 


-0.69 


tpl 


F1315+4438 


2.77 


+ 0.16 
-0.28 


-22.49 


SB2 


<0.08 




0.89 


0.17 


<-1.40 


dist 


F1355+3607 








<0.23 




3.48 


0.17 


<-1.58 


dist 


F1430+3557 


3.5 


+ 0.6 
-0.2 


-24.82 


SB3 


2.13 


0.08 


8.54 


0.17 


-0.81 


unres 


F1445+2702 










13.14 


0.04 


31.07 


0.17 


-0.50 


unres 


F1447+1217 


4.70 


+ 0.13 
-0.18 


-25.11 


SB3 


9.15 


0.14 


48.90 


0.17 


-0.98 


unres 


F1451+0556 


2.567 


+ 0.005 
-0.005 


-26.76 


SB5 


2.57 


0.10 


2.35 


0.17 


0.05 


unres 


F1505+4457N 


1.02 


+ 0.03 
-0.03 


-20.45 


Sa 


1.76 


0.14 


7.84 


0.17 


-0.87 


dbl 


F1505+4457S 


1.02 


+ 0.03 
-0.03 


-20.45 


Sa 


8.04 


0.14 


9.24 


0.17 


-0.08 


dbl 


F1505+4457tot 


1.02 
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dio galaxies able to be used with NIR adaptive optics. 
These 58 objects were observed in B, V, R, I, z', J, H, and 
K bands and their redshifts were estimated using SED 
fitting and generally confirmed as accurate with spec- 
troscopy. 

It was found that the FIRST-BNGS sample objects 
gener ally follow the IR Hubble diagram for radio galax- 



ies. (El Bouchefry fc Cress| ( |2QQ 7) did a study with FIRST 
galaxies in the JNOAO Deep- wide Field Survey Bootes 
field and found similar results. They had a large spread 
at low redshift (z < 0.7) that we did not observe probably 
due to our selection preference for high redshift objects. 
The low redshift sources in our sample would be the 
brightest; we chose not to complete these in all bands in 
favor of fainter objects and as a result we would not have 
photometric redshifts for them. Several objects at high 
redshift (z>l) have best-fit SEDs consistent with young 
stellar populations (SB 1-6). The few long- wavelength 
observations tend to favor either a quiescent or M82-like 
star-forming SED (rather than LIRG or ULIRG type 
SEDs). This may suggest that at least some of these 
high redshift galaxies may be in an active star-forming 
phase which is not what is seen in the K-z Hubble di- 
agram which shows something like a passively evolving 
population from very high redshift (z>5). No trend with 
radio spectral index was found though there were more 
steep-spectrum sources observed than flat-spectrum. 

Many of these objects are at significant redshift, and 
this sample provides a unique tool to study galaxies 
at high redshift. Today's ground-based instrumentation 
even allows sufficient resolution to measure the funda- 
mental parameters of the host galaxies; imaging pro- 
vides resolved morphology and color gradients and spec- 



troscopy allows the dynamics of each galaxy to be stud- 
ied. A recently completed study of the morphologies of 
11 FIRST-BNGS galaxies from a subsample observed 
with the Subaru Telescope suggests that these objects 
tend to be compact, blue, dynamically-relaxed galaxies 
(Stalder & Chambers, in prep.). These intriguing objects 
provide a glimpse of the detailed picture of the first few 
Gyr of the history of the universe that future projects 
such as TMT and JWST will provide. They also hold 
great potential in studying the high redshift universe and 
deserve further attention. 

We thank Michael Connelley, Steve Howell, Elizabeth 
McGrath, and Barry Rothberg for assisting in some of 
the imaging and spectroscopy observations for this huge 
data set. We also acknowledge the telescope support 
staffs at the University of Hawai'i 2.2-meter as well as 
the Infrared Telescope Facility, which is operated by the 
University of Hawai'i under Cooperative Agreement no. 
NCC 5-538 with the National Aeronautics and Space 
Administration, Science Mission Directorate, Planetary 
Astronomy Program. This research was partially sup- 
ported by the Extragalactic and Cosmology division of 
NSF under grant AST 0098349 and also partially sup- 
ported by the Pan- STARRS Camera Group. This re- 
search has made use of the NASA/IPAC Extragalactic 
Database (NED), which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, un- 
der contract to the National Aeronautics and Space Ad- 
ministration. This work is partly based on observations 
made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of 
Technology under a contract with NASA. 



APPENDIX 
SED TEMPLATE LIBRARY 

For objects with m ore than 4 bands of pho tometric data, photometric redshifts were derived using the public SED 
fitting code, Hyperz ( |Bolzonella et al. 2QQQ[ ) which requires a set of templates to fit to the photometric data. See 
section 4 for a detailed descriptions of the overall pro cedure. 



The SED template library used were originally from Kinney et al. (1996), E, SO, Sa, Sb, Sc and SB1 to SB6 (starburst 
models over a range of reddening) b uilt using IUE and opti cal data with a range of wavelength from 1200A to 1.0/im. 
This library was later extended by |Mannucci et aT] ( |2QQ1| ) into the NIR (2.4/i m) using averaged g round-based NIR 
spectra of local prototypical galaxies which closely matched the templates from |Kinney et ah] ( |1996[ ). 

Our high redshift SED fitting procedure requires even shorter wavelength data, and the Spitzer and submillimeter 
data require longer wavelength data than provided by th ese libraries. Therefore f urther UV and mid-IR extensions were 
accomplished using a procedure si milar to that used by Bolzonella et al. (2000) to extend [Coleman et al.| ( [l980| ) tem- 
plates. An appropriate GISSEL98 ( jBruzual fc Charlot|l993p synthetic blue spectrum (constant SFR, and age=0.1Gyr) 
was chosen for Sa, Sb, Sc and SB1-6 or a red spectrum (with delta-function starburst, age=19Gyr) for E and SO. These 
models were chosen to match the overall slope of the continuum where it overlaps with the template. These spectra 
were then grafted onto the Kinney/ Mannucci templates by matching the average continuum levels. The matching 
wavelength ranges were chosen to be relatively smooth (few spectral lines) and flat in F\ (1230 to 1500 A in the UV 
and 2.30 to 2.36/im in the Near-IR). Figures A.l and A. 2 show the full wavelength-range SEDs. 

In order to fit SEDs to Spitzer or submillimeter data, a further extension was applied to these templates using 4 
additional FIR-submm template SEDs chosen to reflect an assortment of possible prototypica l spectra. The templates 
used were Arp 220 ( Bressan et al.|2002 ) representing a ULIRG SEP, M8 2 ( |Bressan et al.|2002|) representing a starburst 
SED, and 2 templates in the synthetic library from Dale et al. (2001) representing LlKGs (a = 1.06) and quiescent 
(a = 2.5) SEDs. The extension was done almost identically to the UV/IR extensions except using the continuum 
levels from 1.95 /im to 2.05 jam to scale the spectra. All 4 templates were applied to each of the SED templates so it 
is possible to fit the SEDs continuously from 100 A to 1 mm. Figure A. 3 shows the 4 FIR-submm extensions to the 
SB3 optical-IR template. 
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FIG. A.l — Template Spectra used in Hyperz from |Mannucci et al~| ( 200 1 ) and extended with GISSEL98 ( jBruzual fe Charlot|1993| ). 
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